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^ . Abstract 
CN . 

' We compare the predictions of different models for transverse momentum spectra in 

■ relativistic heavy ion collisions. Particular emphasis is given to hydrodynamic flow with 

. different assumptions for the radial expansion, and to models based on a superposition 

of fireballs along the transverse direction (random walk). It is shown that present 
experimental data cannot distinguish between the models under consideration. 



1. Introduction 



>< : Recent results on particle yields in Pb - Pb collisions at CERN ^ show that the 

I average transverse momenta of hadronic secondaries are considerably larger than those in 
p — p collisions, and that the increase depends strongly on the mass of the hadron. One 
interpretation of this effect is that the observed hadrons are subject to transverse flow 0; 
this is a velocity effect and therefore heavy particles acquire more momentum than light 
particles. Recently, however, renewed attention has been drawn to the fact that there is a 
"normal" pt-broadening observed in all reactions involving nuclear targets, from Drell-Yan 
dilepton production to low pt mesons or baryons ^]. For high pt hadrons this is generally 
referred to as Cronin effect 0. Here, as well as in Drell-Yan dilepton or in quarkonium 
production, it is accounted for by the fact that successive parton scatterings rotate the 
collision axis relative to the beam axis: any given transverse momentum distribution will 
appear broadened when it is measured in the reference frame fixed by the incident primary 
beams. 

These considerations were applied to lowp^ hadron production in nuclear collisions 0, as- 
suming that successive collisions in nuclear reactions lead to a random walk in the transverse 
rapidity plane. The displacement per collision in transverse rapidity was determined from 
p — A collisions; the normalized pt spectra for A — B collisions are then predicted parameter- 
free and agree quite well with the mentioned data from Pb — Pb interactions lH, 0. In 
particular, this "normal" pt-broadening also reproduces the increase with increasing hadron 



mass, giving more broadening for nucleons than for kaons, and more for kaons than for pions. 
A very recent study 0] has gone even further and determined the "kick per colhsion" from 
p — p rather than from p — A data. 

The Pi-broadening observed in nuclear reactions can thus be quite well accounted for 
by random walk collision axis rotations. Perhaps one might consider such a phenomenon 
as a precursor for "transverse flow". Nevertheless, any hydrodynamic description of pt- 
distributions from A — B collisions, with the flow velocity as open parameter, has to face 
two questions: why is there also broadening in p — A interactions? and why can the "flow 
velocity" in a random walk approach be determined from p — Aoi even from p — p collisions? 
Perhaps only two-particle correlations, rather than single particle spectra, can distinguish 
between hydrodynamic flow and a random walk approach |p. 

In this paper we want to show explicitly that the present information on transverse 
momentum spectra can be adequately explained in conceptually different models and thus 
cannot determine the origin of the observed broadening. 

We will first consider a hydrodynamic model with transverse flow and rapidity scaling. In 
this case, all spectra show a characteristic flattening for very small values of pt- This is due 
to the flow which adds momentum to very slow particles and hence depletes this kinematic 
region. Next, we will study two simplifled forms for the transverse velocity: a constant 
velocity and one with a polynomial dependence on the transverse radius. All three versions 
are able to reproduce the data, with some caveats to be discussed later on in the paper. In 
particular, the degree of flattening characteristic of a full hydrodynamic picture does not seem 
to be present in the data. As alternative description, we consider the behaviour obtained 
from the random walk model of . In this case, one typically overshoots the data somewhat 
at low transverse momenta. This could well be due to the gaussian distribution used to 
describe the successive kicks of the projectile; it has the tendency to accumulate particles 
around the origin, and an exponential distribution might lead to a better description of the 
present data. Once the quality of the data improves, such a flt should deflnitely be made; in 
the present paper, we restrict ourselves to the models so far proposed in the literature. An 
inherent difficulty of the random walk model is that each one of the successive collisions is 
determined by results obtained from p — A collisions, so that in particular the relative particle 
abundances will also be the same asvnp — A collisions. There is therefore no mechanism to 
obtain the observed increase in the abundance of strange particles. 

2. Hydrodynamic Expansion 

The momentum distribution of particles is given by the well-known Cooper-Frye for- 
mula 



where the integration has to be performed over the freeze-out surface described by cr^. For 
the temperatures under consideration it is safe to neglect quantum statistics and we will 
therefore work with the Boltzmann distribution from now on; the generalization to Fermi- 
Dirac or Bose-Einstein statistics is straight-forward. We thus have 




(1) 



f{x,p) = exp[(-p.M + /i)/T], 



(2) 
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where T, /z and are the (space-time dependent) temperature, chemical potential and 
four-velocity of the fireball, respectively. For a fireball at rest we have 

w'^ = (l,0), (3) 

for one boosted in, e.g., the x-direction, 

■u'^ = (cosh?/r, sinh?/T, 0, 0), (4) 

while for a boost in the azimuthal direction cf) we get 

= (cosh ut, cos sinh ut, sin (j) sinh yx, 0). (5) 

For a static fireball undergoing instantaneous freeze-out, the corresponding freeze-out surface 
is given by 

da^" = {d^x.Q), (6) 
so that we obtain the standard expression (in Boltzmann approximation) 

For boost-invariant cylindrical expansion along the z-axis (recall d'^x = TdTdyrdrd(f)), we 
have 

da^ = {Tdyrdrd(f),TdTdyrd(f)), (8) 

where the second component is in the f direction, i.e. perpendicular to the surface of the 
cylinder. For the case where the flow is azimuthally symmetric, i.e. when an average is made 



over all events or when only head-on collisions are considered, one has therefore [jTO|, p 
f dN \ 
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where Tpir^ refers to the freeze-out time which in general depends on r, so that the center 
of the cylinder freezes out before the surface. 

We have taken the values for the derivative dTp{r)jdr from reference [|T3|; in the figures, 
the results based on eq. (9) are labeled ''hydro 2+1 dim". 

3. Transverse Flow 

If the space-time development typical of the hydrodynamic expansion is not taken into 
account and one instead considers bubbles of fluid receiving boosts in the transverse direction, 
then it is natural to take the freeze-out time as independent of r, 

Tpir) = Tp (10) 



so that everywhere within the volume the particles freeze out simultaneously and the second 
term in equation (9) disappears. One is then left with 
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where V is defined to be nRpTp, where Rp denotes the value of the radius at freeze-out. 
The variable ^ = r/Rp. is directly related to r, measuring the distance from the axis of the 
cylinder. A natural possibility is to allow for a velocity profile of the form 
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tanh|/4 = vt 



,a + 2 



r 
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(12) 



In the figures we will label the results based on equation (11) (with a = 1) as "f ^-profile" . 

If the transverse expansion is taken to be independent of r, so that the expansion velocity 
is the same in the center as it is on the surface, then one obtains 
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(13) 



where the volume V is again determined hy V = ixR^Tp. It should be noted that a constant 
transverse velocity leads to a problem at the origin, since the transverse velocity must be 
zero at this point because of symmetry considerations. In the figures we will refer to the 
results of eq. (13) as "vt=const" . 

4. Random Walk 

In the random walk approach, nuclear collisions are assumed to be much like elementary 
p — p collisions, except that in the successive scatterings occurring in nuclear targets, the 
collision axis will be rotated. 

It was shown recently |]I4| that a thermal description accounts quite well for the particle 



ratios in — e~ , p — p and p — p collisions. We therefore follow this picture and assume 
that in each of the successive interactions of a nuclear collision one creates a fireball just like 
that formed in a nucleon-nucleon collision; there is no need to introduce hydrodynamic flow 
for such little fireballs. The only difference now is that after the first collision the next one 
will generally occur at some non-vanishing transverse velocity. We therefore have to know 
the propagation of transverse momentum through successive collisions. It seems simplest to 
assume that this will follow a random walk pattern, which in |^ was taken to be Gaussian. 
As mentioned, one should eventually study different distributions in such an approach; we 
will restrict ourselves to the Gaussian 



fpAip) = 

where p is the transverse rapidity and 
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exp(-p /6p^), 



(14) 



= {Na - 



(15) 
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denotes the kick per collision S as determined from p — A interactions. The corresponding 
distribution for an A — B collision can then be obtained by the convolution 



fAB{p)= / fpA{p')fpB{p")S{p' + p" - p)dp'dp", 



leading to 



with 



fAB{p) 
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exp(-pV'^Ai?) 



(16) 



(17) 



^Ib = {Na + Nb- 2)5\ (18) 

For simplicity we have followed the analysis of and have taken all fireballs at the same tem- 
perature; again this can easily be generalized to a more general situation with a distribution 
in temperature. The final expression thus becomes here 
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It should be noted that the volume in the eq. (19) refers to the volume of the system as 
observed in a p — p collision, since each collision in the random walk produces a. p — p type 
of fireball. For normalized distributions, it will of course drop out. 

If we now introduce a boost-invariant distribution of fireballs along the longitudinal rapi- 
dity axis, we finally obtain by integrating over the fireball distributions 
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In the figures the results based on eq. (20) will be denoted by "random walk' 
5. Results 



(20) 



We now want to compare the models described in the previous sections with the latest 
experimental data from Ph — Pb collisions obtained by the NA44 and the NA49 [0 
collaborations at CERN. These data were obtained in different kinematic regions. In the 
figures, we have simply scaled the NA44 data to the NA49 data. The parameters used in 
the different descriptions are listed in Table 1. 

In Fig. 1, we show the NA49 data for negative hadrons together with the NA44 results 
for 7r"^. All models show good agreement with the NA49 data, but they cannot reproduce 
the steeper NA44 vr^ behaviour at small pt together with that in the larger NA49 pt range. 

In Fig. 2, we show the NA49 data for the surplus of positive hadrons, — h~ , and the 
NA44 data for protons. In this case there is good agreement between both data sets and all 
models, except that random walk overshoots and the full hydrodynamic model undershoots 
the behaviour at very low pt. As noted, in the random walk picture, the choice of another 
distribution in the "kick" -parameter 6 may well bring the spectrum down. The characteristic 



A 



low pt dip in the hydrodynamical model is reduced somewhat when we choose a different 
parameterization of the transverse velocity. 

In Fig. 3, we show the NA49 data for neutral hyperons, i.e. A and In this case, the 
hydrodynamical description with constant Vt does quite poorly, while all other models give 
a reasonable description of the data. 

In Fig. 4, we show the NA49 data for K'^ and the NA44 results for K~^. There is 
again good agreement between data and models, except in the low pt region, where the full 
hydrodynamic model shows a depletion not present in the measured yield. 



Model 


T/ [GeV] 


< f _L > 


References 


Vt = const 
Vt — profile 
random walk 
hydro 2+1 dim 
hydrodynamics ([|19]) 
hydrodynamics ([pO|) 
from particle ratio 


0.15 

0.12 

0.15 
0.115 

0.14 
0.143 
0.17 - 0.19 


0.37 
0.43 

0.5 
0.34 


[12] 
[6] 

|iq 

m 

[14,15,17,18] 



Table 1: Results for freeze-out parameters in different models. 

Finally we show in Fig. 5 how the different fits tend to disperse for larger values of the 
transverse momenta. A similar behavior is obtained for all particle spectra. 

6. Conclusions 

We have compared the latest experimental results on transverse momentum spectra in 
Pb — Pb collisions with several different models. All of them manage to describe these data 
reasonably well, apart from some low pt discrepancies; hence transverse spectra so far really 
do not allow us to decide between a hydrodynamic or random walk origin of the observed 
broadening. 

Moreover, several further problems remain unanswered: all models call for a rather low 
freeze-out temperature, as seen in Table 1, where we have for comparison also listed results 
obtained in refs. [|l^, Our results are typically in the range Tp ~ 120 — 150 MeV. This 
is low when compared to the temperatures needed to explain the hadron abundances in 
particle ratios, where one needs values in the range of 170-190 MeV but it must be noted 
of course that the various hydrodynamic models do not take into account resonance decays. 
The random walk model describes the data quite well but it will lead to the same hadronic 
abundances as those observed in p — A collisions. Thus it cannot reproduce the increase in 
strangeness production reported by all experimental groups. Besides the uncertainty about 
the origin of the observed transverse broadening, we thus also still lack a consistent account 
of both slopes and particle abundances. 
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Figure Captions 

Figure 1 : Comparison between model predictions and the data for pion distributions. 
Figure 2 : Comparison between model predictions and data for proton distributions. 
Figure 3 : Comparison between model predictions and data for the A distribution. 
Figure 4 : Comparison between model predictions and data for kaon distributions. 
Figure 5 : Comparison between model predictions and data for the deuteron distribution. 
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